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(N : ABSTRACT 

' The tight correlation between the masses of central black holes and their host 

^ [ spheroids in nearby galaxies and active galactic nuclei (AGN) suggests that black 

in 

(N 



hole growth is closely related to their spheroid formation. Based on our previous 
work regarding such a joint evolutionary scheme and the consequential black hole 
I to bulge mass correlation, we use the X-ray luminosity function of AGN and the 

cosmological evolution rate which are from ROSAT X-ray surveys to estimate 



O 



Qh| the cosmic star formation history associated with the black hole growth. By the 

O ' basic assumption that the major black hole growth occurs during the luminous 

^ ' AGN phase, the luminosity function of AGNs as a function of redshift traces not 

' only the accretion history of the black holes but also the cosmic star formation 

history of the spheroids. 

^ ' Although the space density of the especially luminous QSOs is very low, we 

■ show that the total amount of star formation associated with the massive black 

hole growth is almost the same as that of Lyman Break Galaxies detected by the 
current optical deep surveys. We thus argue that the optical deep surveys may 
miss about half of the net star formation in our Universe. This is probably due 
to in part significant dust extinction as well as the small field of view of previous 
optical surveys which cannot sample such rare events with relatively short time 
scale. However, the far infrared emission from the dust heated by star formation 
on-going during the black hole growth could sufficiently account for the observed 
SCUBA number counts, and would be the probable dominating energy source of 
the SCUBA population. 



Subject headings: galaxies: elliptical and lenticular, cD-galaxies: active-galaxies: 
evolution-quasars: general-(ISM:) dust, extinction 
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1. Introduction 

Recent observations reveal a tight correlation between the masses of central black 
holes and their host spheroids in nearby galaxies and AGNs (Kormendy & Richstone 1995; 
Faber et al. 1997; Magorrian et al. 1998; Laor 1998; Ho 1999; Wandel 1999a). A mean 
value of the black hole to bulge mass ratio, Mbh/Mbuige ~ 10~^ and the local black hole mass 
density, pbh ~ 3 — 5 x 10^ Mq/Mpc^ are now fairly agreed upon by observations using differ- 
ent techniques (Salucci et al. 1999; Kormendy & Gebhardt 2001; Merritt & Ferrarese 2001; 
McLure & Dunlop 2002). This correlation implies a possible scenario where the black hole 
accretion history in the inner few parsecs and the ~ kpc scale star formation during spheroid 
formation are closely connected. Simply speaking, the amount of gas accreted to grow the 
central black hole is proportional to the amount of gas driving the star formation to popu- 
late the spheroid within a similar time scale of the two activities (Wang & Biermann 1998; 
Wang et al. 2000; Haehnelt & Kauffmann 2000; Monaco et al. 2000; Burkert & Silk 2001; 
Page et al. 2001). 

This relation also opens an important avenue to the understanding of the cosmic star for- 
mation history (CSFH) associated with the black hole growth by accretion, complementary 
to the previous optical/near- infrared (NIR) or far- infrared deep surveys. 

The study of the CSFH from the optical/NIR deep surveys describes only the un- 
obscured star formation in galaxies, and may be blind to activities during the massive 
spheroid formation because of the significant dust extinction or the small field of view 
(Madau et al. 1996; Steidel et al. 1996; Connolly et al. 1997; Steidel et al. 1999). In fact, 
our knowledge of the cosmic star formation history, especially at high redshift, has dra- 
matically changed in the last few years due to the results of the far-infrared and sub- 
millimeter deep surveys. The compelling lines of evidence show that the absorption and 
re-radiation of light by dust in the early epoch of galaxy formation and evolution is sig- 
nificant (Puget et al. 1996; Guiderdoni et al. 1997; Schlegel et al. 1998; Fixsen et al. 1998; 
Hauser et al. 1998; Sanders 1999; Barger et al. 1999). 

Although the observations at far-infrared and submillimeter wavelengths may help to 
unveil the early dusty star formation epoch, it is still a very difficult job to obtain an unbiased 
view of the early star formation history from these luminous infrared systems due to the 
rather poor angular resolution of SCUBA and to the faintness of the optical counterparts 
(Blain et al. 1999a; De Zotti et al. 2001). 

The main purpose of this paper is to fairly reconstruct the global star formation history 
associated with AGN accretion and give a proper estimation of the star formation activity 
during spheroid formation. We assume the X-ray luminosity Lx of the AGN is powered by 
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accretion onto a central massive black hole, especially during the luminous phase. Deep X- 
ray surveys provide a direct probe of the AGN accretion history, and the history of the black 
hole growth as well as the joint star formation, based on the assumption that the luminous 
AGNs reflect the stage of major black hole growth and spheroid formation (Boyle et al. 2000; 
Miyaji et al. 2000). In this paper, we trace the AGN evolution with the soft X-ray 0.5 — 
2 keV local luminosity function and the evolution rate by Miyaji et al. (2000), which gives 
an excellent fit to the ROSAT surveys of various depth, including the number counts and 
the soft X-ray background. 

The X-ray view of the CSFH is recently discussed by several authors for normal spi- 
ral galaxies where the X-ray emission is dominated by a population of X-ray binaries, 
hot interstellar gas, or even low luminosity AGNs (LLAGNs) with fiuxes about 10~^^ — 
10-^^ erg cm-2 g-^ at 0.5 - 2 keV energy band (Cavahere et al. 2000; Ptak et al. 2001; 
Hornschemeier et al. 2002; Miyaji et al. 2002). It is not the purpose of this paper to give 
any constraints on the star formation in galaxies or LLAGNs. Instead, we study the inten- 
sive star formation associated with black hole growth in massive spheroids using the X-ray 
surveys. We adopt in the calculation: 1) a simple connection between the X-ray emission 
from AGNs and the black hole mass by assuming an Eddington ratio e; 2) the black hole to 
bulge mass correlation; 3) a similar time scale for the black hole growth and the intensive 
star formation populating finally the spheroids based on our previous work. The present 
day black hole density and the SCUBA number counts are used as two important model 
constraints in the calculation, with the set of cosmological parameters (Qm, Q\) — (0.3, 0.7) 
and Ho = 50 km/s/Mpc. 

2. Black hole mass distribution and the accretion history from ROSAT 

surveys 

The mass accretion onto the central black hole is an efficient source of X-ray radiation, 
especially for the luminous AGNs. Although X-ray binaries can contribute as much as 
~ 5% - 10% of the X-ray background (XRB) fiux in the 0.5 - 2 keV band, the bulk 
of the energy density of the XRB is certainly explained by AGNs (Mushotzky et al. 2000; 
Barger et al. 2001a; Giacconi et al. 2001). Thus, the deep X-ray surveys could present good 
opportunities for the study of accretion history of AGNs. We estimate the black hole mass 
from the X-ray luminosity by assuming an Eddington ratio e, which is defined as the fraction 
of Eddington luminosity at which AGN radiate (e = ^^^)- With the bolometric correction 
/3 (Lboi = /3Lx), we find the black hole mass Mbh = olfc- Lx is the AGN 0.5 - 2 keV 
luminosity in units of lO^'^ erg s~^ and Mbh in units of Mq. The black hole (BH) mass 
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function at different redshifts is converted from the X-ray luminosity function $(i^x, z) by 
assuming reasonable values for the Eddington ratio e and the "duty cycle fon" of the AGN 
active phase. Following the findings of various AGN observations, e depends weakly on 
the luminosity, where most luminous QSOs radiate at about the Eddington limit and low 
luminosity AGNs (L < 10^^ erg/s) show e ~ 0.1 — 0.05. We thus approximate e = 
lQ^(iogL-49). ^ ^ Q 2 is a scaling factor (Padovani 1989; Salucci et al. 1999; Wandel 1999). 
The "duty cycle fon", i.e. the fraction of black holes which are active at a given time of 
redshift z, may be a function of several parameters, such as the redshift and the quasar 
light curves. In this work, we simply assume that each AGN shines for a constant time tg, 
fon = tQ/tHub(z), where tHub(z) is the redshift dependent Hubble time. This means that we 
simplify the QSO activity as a single significant burst with shining ratio e and duration tq 
(Haiman & Menou 2000). The QSO life time tq = 5 x 10® yrs is adopted in this calculation, 
which is close to the e-folding time te = Mbh/M = 4 x lO^r^/e shown in theoretical models 
with the radiation efficiency t] ~ 0.1, and consistent with the new results of Chandra for the 
QSO accretion duration (Wang & Biermann 1998; Burkert & Silk 2001; Barger et al. 2001). 

Although there are a variety of AGN spectra in the current sample of ROSAT X-ray 
regime by Miyaji et al. (2000), non type 1 AGNs are only a small fraction of the total sample 
and excluding them does not change the main results significantly. In this case, we consider 
all AGNs in the sample of ROSAT soft X-ray regime by Miyaji et al. (2000) as unobscured 
type 1, and include a bolometric correction /3 to convert X-ray luminosity Lx to black hole 
mass Mbh- P — 20 is adopted in our calculation based on the mean type 1 AGN spectral 
energy distribution from Elvis et al. (1994). The black hole mass function ^^^^j^^'^^ can be 
derived from the observed X-ray luminosity function ^^l^"^^ 



A Luminosity Dependent Density Evolution(LDDE) model is suggested by Miyaji et 
al.(2000), where the evolution rate drops with decreasing AGN luminosity, and provides a 
good representation of the available X-ray observational constraints. We follow the LDDE 



black hole mass density accreted during the active phase with lookback time. 

However, there are several lines of arguments which suggest that the optical and soft 
X-ray surveys may miss a large number of type 2 AGNs, especially the high luminosity 
type 2 QSOs at high redshift (Lx > 10^^ erg s~^). Synthesis models of the X-ray back- 
ground in particular require a large number of luminous obscured objects to reproduce the 
2 — lOkeV source counts at relatively bright fluxes (~ lO"^'^ erg cm~^ s~^). The ratio of 



d $(^, A4h) _ 0.013 e d ^{z,L^ 
d Mbh ~ Tf^n dLx 



(1) 



cosmological evolution for the ROSAT soft X-ray luminosity function 



and trace the 
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type 2 to unobscured type 1 Seyferts (Lx < 10^ erg s'^) and the luminous type 2 to 
type 1 QSOs are usually in the range of 4 — 10, in order to get a best fit of the avail- 
able X-ray observations (Maiolino & Rieke 1995; Fabian & Iwasawa 1999; Gilli et al. 2001; 
Maiolino et al. 2001). Although the Chandra and XMM- Newton deep surveys have recently 
detected several examples of the type 2 QSOs, Alexander et al.(2001) suggests that their 
abundance is definitely much lower than what is expected by the XRB population synthesis 
models, and probably in a ratio of about 1 ~ 2 (Norman et al. 2002; Stern et al. 2002; 
Hasinger et al. 2001; Akiyama et al. 2002). 

Since the Miyaji et al. (2000) sample may represent only the population and evolution 
of the type 1 AGNs (i.e. type 1 Seyferts and type 1 QSOs), we adopt a reasonable fraction 
of obscured type 2 sources in our calculation to include the amount of star formation in the 
host galaxies of the obscured objects. Maiolino & Rieke (1995) derived an estimate for the 
ratio of the local absorbed to unabsorbed low luminosity AGNs around a value of 4. Unlike 
the absorption distribution of local Seyferts, the existence and abundance of the luminous 
type 2 QSOs is very uncertain. The difficulties are that: l)some type 2 QSOs could be 
hidden in the Ultraluminous Infrared Galaxies (ULIGs) (Kim & Sanders 1998); 2)although 
it is still to be verified, some type 2 QSOs appear as normal blue, broad-lined QSOs in the 
optical (Halpern et al. 1999). 

However, the present observational data suggests that the evolution of the obscured 
objects might be different from that of the unobscured ones, with the number ratio of 
type 2 to type 1 QSOs increasing with rcdshift (Gilli et al. 1999; Pompilio et al. 2000; 
Reeves & Turner 2000; Gilli et al. 2001). In this case, we divide the soft X-ray Luminosity 
Function(XLF) into two luminosity regions according to the 0.5 — 2 keV e-folding luminosity 
Lg = 10^^-^ erg s~^. The abundance ratio of type 2 to type 1 Seyferts (Lx < Lg) is set equal 
to 4 and the ratio of type 2 to type 1 QSOs is simplified as a (1 -|- zy. a and p are free 
parameters. We also explored the possibility that the abundance ratio of LLAGNs is much 
higher than 4, i.e. as high as 10. It seems that the current far-infrared and submillimeter 
deep surveys are not robust enough to contrain this value. Probably, it is because LLAGNs 
do not contribute much to the bright SCUBA number counts. In this case, the abundance 
ratio of the low luminosity Seyferts is not critical for this work. Thus, the abundance ratio 
of type 2 to type 1 AGNs including low luminosity Seyferts and QSOs can be described as 
the following: 

i?2_i = 4e-^ + a(l + ^)*'(l-e"fe) (2) 

Considering the present day black hole density pbh ~ 3 — 5 x 10^ Mg/Mpc^ and the 
submillimeter deep surveys as the model constraints, we found the ratio of the obscured 
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type 2 to unobscured type 1 QSOs has an upper limit of about 2 {a and p have the best fit 
values of 0.1 and 1.9). There is not much room left for a ratio far beyond this number. This 
is actually consistent with the results of recent Chandra deep surveys by Alexander et al. 
(2001) who suggests that this ratio is unlikely to be more than 8 and is probably considerably 
lower. Since X-ray emission directly measures the accretion of AGNs, we show in Fig. 1 the 
accretion history from the model calculation. The solid line illustrates the accretion history 
of the type 1 AGNs only; while the dashed hne is for the case with type 2 AGNs included. 
The accretion rate density pbh in units of Mq yr~^ Mpc""^ from a multiwavelength study of 
69 hard X-ray selected sources in the field of Abell 370 cluster (A370) , Hawaii Survey Field 
(SSA13) and Chandra Deep Field-North (CDF-N) for four redshift bins are also ploted in 
Fig. 1, where the blue asterisks represent the accretion onto the AGNs calculated from their 
bolo metric luminosities and the red crosses are from their X-ray luminosities as a low limit. 
Fig. 1 shows that the predicted accretion rate density of the type 1 AGNs is below the upper 
bound derived from the hard X-ray selected samples (Barger et al. 2001). Meanwhile, we see 
that the fraction of type 2 AGNs included in our model still agrees with the current results 
of deep hard X-ray surveys. It is worth mentioning that the treatment of the abundance 
of type 2 AGNs in this work is very rough. We need more data from the Chandra and 
XMM-Newton surveys at significant fainter fluxes to better constrain the abundances. 

3. Star formation associated with black hole growth and the SCUBA counts 

Recent observations seem to support an evolutionary scenario where the black hole 
growth and intensive star formation during the spheroid formation are co-evolving events. 
In this scheme, the accretion onto a seed black hole and the rapid star formation in a ~ kpc 
region are significantly enhanced due to a merging process or tidal interactions. Accretion 
-|- star formation compete for the gas supply and may self-regulate by the fundamental 
properties of the viscous accretion disk, leading to the formation of massive black holes 
scaling with their spheroids. In this case, the mass function of spheroids may have a similar 



form to the balck hole mass function, and the spheroidal mass distribution 
be derived by: 




d ^z, M,ph) ^ d'l'(z,Mbh) 
d Mgph d Mbh 



. R(Mbh/M,ph) 



(3) 



R(Mbh/Msph) is the black hole to bulge mass ratio. We adopt here a mean value R(Mbh/Msph) ~ 

0.002 as the first approximation (McLcod & Riekc 1994; Kormendy & Richstone 1995; Faber et al. 1997; 

Magorrian et al. 1998; Merritt & Ferrarese 2001a). 
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The aim of this section is mainly to estimate the star formation rate during spheroid 
formation in individual galaxies, and then discuss the co-moving star formation history 
related with such an epoch. The SCUBA number counts are used as a model constraint or 
consistency check for the amount of star formation, especially at z > 1. 

The basic assumption of this calculation is that the stellar population of the spheroids 
arc mostly from an intensive star formation phase in the central kpc region which is asso- 
ciated with the mass accretion onto the central black holes. In the model calculation of 
Wang & Biermann (1998), star formation would dominate the energy output in the first 
~ 5 X 10® yrs for the stellar black hole growing to a mass of ~ 10^ M0, even with full- 
Eddington accretion. Because of the heavy dust extinction, the intensive star formation 
during this process would emit a large fraction of its energy at infrared wavelength. The 
Ultraluminous Infrared Galaxies (ULIGs) may represent such a dust-enshrouded stage of 
both intensive starbursts and AGNs. Once the AGN activity becomes powerful enough, the 
nucleus sheds its obscuring dust, allowing the AGN to visually dominate, completing the evo- 
lutionary process from ULIGs to optically selected QSOs (Sanders et al. 1988; Yamada 1994; 
Mihos & Hernquist 1996; Silk & Rees 1998; Taniguchi et al. 1999). 

We adopt a characteristic star formation time scale Tgf for the starburst activities coupled 
with the black hole growth, which is approximately equal to the QSO life time tq (i.e. 
Tsf ~ ~ 5 X 10® yrs). This is consistent with the predictions of the theoretical models 
and the recent Chandra results of the QSO accretion duration (Wang & Biermann 1998; 
Burkert & Silk 2001; Bargcr et al. 2001). The star formation rate during such a phase would 
be SFR = Msph/Tsf. Following the cosmological evolution inferred from X-ray deep surveys, 
we get the redshift dependent black hole mass distribution by eq. 1, and the mass function of 
spheroids with redshift by eq. 3. So far, we get to know the star formation rate in individual 
spheroid during the past active phase, the mass distribution and the cosmological evolution. 
The co-moving star formation rate vs. redshift from our calculation is shown in Fig. 2. 

In section 2, we have applied the local black hole density to constrain the black hole 
mass distribution and cosmological evolution in our calculation. Independently, submillime- 
ter deep survey would be a good tracer and constraint of such an active star formation 
stage associated with the major phase of black hole growth and spheroid formation. To 
discuss the contribution of such an intensive star formation phase to the SCUBA number 
counts, we should convert the star formation rate to far-infrared emission for individual 
source. This relation is not straightforward even for the local ULIGs because of the high 
internal extinction. A number of authors have discussed how the star formation rate in a 
galaxy can be inferred from its optical, UV or far-infrared luminosity (Scoville & Young 1983; 
Leitherer & Heckman 1995; Rowan-Robinson et al. 1997). The transformation factor de- 
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pending upon star formation parameters could vary by more than one order of magnitude. 

Wc adopt a mean value Lfir/Lq = 3.8 x 10^ SFR (MqYI"^), and a mean color ratio 
Rc = Lfir/L85o ~ 5 X 10^ given by Chary & Elbaz (2001) from IRAS, ISO and SCUBA 
surveys. 

Similar to eq. 1, we can derive the 850 //m luminosity function at different redshift as 
follow: 

d$(^,L85o) ^ d<!>{z,M,p^) dM,ph 

d 1/850 d -^sph d -^850 

where f^i = tsf/tnub reflects the fraction of galaxies which are in the active stage with 
intensive star formation ongoing. = 3T>rw^L^' and tgf are the parameters adopted 

in our model and discussed above, u — 3.5 x 10^^ is the frequency of 850 fj,m. Since 
early star formation could be well constrained by the far-infrared and submillimeter deep 
surveys, we show the 850 number count fitting in Fig. 3. We find that the amount of 
star formation related to the spheroid formation during the black hole growth heats the dust 
sufficiently to account for the far-infrared emission in most of the SCUBA counts. 

Throughout the discussion above, we only focus on the star formation history in massive 
spheroids, i.e. elliptical galaxies or the disk galaxies with bulges. Although there are two 
kinds of bulges in the disk galaxies strengthened by especially Hubble Space Telescope ("mini- 
elliptical" vs. " pseudobulge" ) and their formation mechanism might be very different, both 
of them follow the same correlation of Mbh — MB,buige or Mbh — Oe, presumably consistent 
with a scenario that the bulge formation and the black hole growth are closely connected 
(CaroUo et al. 1997; CaroUo et al. 1998a; CaroUo et al. 1998b; Peletier et al. 2000). The 
star formation during the massive spheroid formation would be drastic and dusty, very 
probably missing in the current UV/optical deep surveys. However, some disk galaxies 
harbor only small black holes with the black hole to bulge mass ratio much lower than that of 
the spheroidal systems. The typical example is M33, a bulgeless galaxy with the upper limit 
on a black hole mass Mbh < 1000 Mq (three magnitude lower than what is expected from the 
mass correlation in massive systems). The star formation in these galaxies may be mild and 
the dust extinction is not so severe. Normally they are contained in the optical deep surveys, 
or other studies of the spiral galaxies (Cavaliere et al. 2000; Hornschemeier et al. 2002). The 
detailed discussion of the star formation history of these galaxies are beyond the scope of 
this paper. 



(4) 
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4. Discussion and summEiry 

We have begun to probe the cosmic star formation history associated with AGN ac- 
cretion by X-ray deep surveys based on the tight correlation of Mbh — Me.buigc in early 
type galaxies and nearby AGNs. This approach is parallel and complementary to the cur- 
rent study from optical and infrared observations of the star formation activities at z > 1 
(Madau et al. 1998; Pettini et al. 1998; Hughes et al. 1998; Barger et al. 1998; Ellis 1998). 

The cosmic star formation history associated with AGN accretion derived from our cal- 
culation is approximately comparable to that of the normal galaxies in Madau plot after a 
reasonable dust correction, although they are much rarer objects compared with optically se- 
lected galaxies (The co-moving number density of luminous AGNs is as small as 10~^ Mpc~^, 
see Fig. 2). In this case, we might say that about half of the star formation (if not more) 
in our Universe is closely connected with AGN accretion, and we might severely overlook 
the intensive star formation during the epoch of the spheroidal formation by optical/NIR 
surveys due to the significant dust extinction and the small sample volume of these surveys. 

We found from Fig. 2 that the peak of the intensive star formation representing the 
spheroid formation is at z ~ 2, not necessarily much beyond this epoch even if we take 
into account a reasonable fraction of type 2 QSOs at high redshift. According to various 
observation, we divide the soft X-ray luminosity function into two luminosity regions, where 
the abundance ratio of the type 2 to type 1 Seyferts is set equal to 4, and the ratio of the 
type 2 to type 1 QSOs is simplified as a power law function of a (1 -|- zy. Actually an 
upper limit of the abundance ratio of the type 2 to type 1 QSOs is about 1 ~ 2 from our 
calculation with the constraints of the local black hole mass density and the results of the 
submillimeter deep surveys. However, current model constraints are not sufficient to reject 
a much higher abundance ratio for low luminosity type 2 to type 1 AGNs (> 4). The reason 
may be that the bright SCUBA counts are dominated by the luminous AGNs and the ratio 
of the low luminosity AGNs is not critical. Future results of the Chandra and XMM/Newton 
deep surveys would give more information on the abundance of the obscured objects. 

The energy budget of the submillimeter sources has been discussed by several au- 
thors, which show that the AGN powered far-infrared emission in the obscured objects 
could account for only a certain fraction of the SCUBA number counts (Almaini et al. 1999; 
Risaliti et al. 2002). The SCUBA/X-ray anti-correlation of the Chandra deep surveys gives a 
clear detection of AGN activities in those samples only ~ 10% (Bautz et al. 2000; Fabian et al. 2000; 
Barger et al. 2001b; Hornschemeier et al. 2001). Nevertheless, Fig. 3 shows that the far- 
infrared emission from the dust heated by the intensive star formation during the black 
hole growth could sufficiently interpret the number counts of the submillimeter deep sur- 
veys, where the SCUBA number counts are dominated by the contribution from the star 
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formation in the host galaxies of the luminous type 1 or type 2 QSOs, and the the faint 
parts (S850 < 1 mJy) may be from those low luminosity X-ray sources with small bulges. 
In this case, our calculation may indicate that the star formation activity might dominate 
the energy power (at least comparable to the AGNs) in the far-infrared emission in these 
SCUBA sources, consistent with the multiwavelength observations of the submillimeter se- 
lected galaxies. They suggest that even when an AGN is present in a SCUBA source, it 
rarely dominates the engergy budget of the galaxy (Frayer et al. 1998; Alexander et al. 2002; 
Ivison et al. 2002; Small et al. 2002). 
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Fig. 1. — Accretion rate density (ARD) in units of Mq yr~^ Mpc~^ vs. redshift from the 
model calculation. The solid line illustrates the calculated accretion history of type 1 AGNs, 
while the dashed line is for the case with a reasonable fraction of type 2 AGNs included. 
The data points, blue asterisks represent the accretion density calculated from the bolometric 
luminosities of 69 hard X-ray selected sources; while the red crosses are calculated from their 
X-ray luminosities as a low limit (Barger et al. 2001). 
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Fig. 2. — Co- moving star formation rate vs. redshift predicted by the model. The data 
sources are indicated in the Figure, while the filled symbols are those with dust correction 
(x4.7) adopted by Steidel et al. 1999. The upper line (green) is the star formation in the 
host galaxies of whole AGN sample including type 1 and a reasonable fraction of type 2 
AGNs, while the lower line (red) is the contribution from type 1 AGNs only. 
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Fig. 3. — 850 fim number count fitting from tfie model calculation in case of type 1 AGNs 
only and the case of type 2 AGNs included. The data are from Barger et al. (1999) and 
Blain et al. (1999). We see from the figure that the star formation powered far-infrared 
emission from X-ray luminous AGNs {L > 10'^^ erg/ s) would dominate the SCUBA counts, 
while the faint parts (Sgso < 1 mJy) might be from those low luminosity X-ray sources. 



